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1 Introduction 



Quite recently, we introduced the minimal N = 2 supersymmetric extension of the Toda chain 
called f-Toda |T[. We have shown that the N = 2 supersymmetric Nonlinear Schrodinger (NLS) 
hierarchy in the N = 2 superspace possesses a discrete group of integrable mappings || which 
is equivalent to the f-Toda chain. Its origin was analyzed in H. In contrast to other chains, 
the f-Toda chain is not algebraically solvable. Such type of chains has not been considered in 
the literature before. 

The goal of the present Letter is to prove the integrability of the f-Toda chain and to 
construct its general solution under appropriate boundary conditions, including the condition 
corresponding to fixed ends. This last problem is equivalent to the problem of constructing 
multi-soliton solutions for each integrable system belonging to the N = 2 super-NLS hierarchy. 



2 The N = 2 super-NLS and f-Toda superfield equations 

In this section, following [lj, we briefly describe the f-Toda chain equations and their relation 
to the N = 2 super-NLS hierarchy. 

Let us proceed with the N = 2 super-NLS equation [|] 

^- = f" + 2D(fjDf), { 5- = -J" + 2D(fjDj) } (2.1) 

where f(x,8,8) and f(x,9,9) are chiral and antichiral fermionic superfields, 

Df = 0, Dj = 0, (2.2) 

respectively, in the N = 2 superspace with one bosonic z and two fermionic 9, 9 coordinates^; 
D and D are the N = 2 supersymmetric fermionic covariant derivatives 

*-!-5£ D 2 =~D 2 = 0, ^ (2.3, 



The chirality condition (|2.2f) reduces the number of independent components of the fermionic 
superfield / (/) from four to two. 

Equations fl2.1|) , as well as each system of the equations belonging to the N = 2 super-NLS 



hierarchy, are invariant with respect to the following discrete f-Toda transformation [|1|]: 

?7- fJ=(HD7Dj))'. (2.4) 
This relation is the definition of mapping, i.e., the rule determining the correspondence between 
two initial superfields / and /, and two final ones, / and /. The notation / (/) means that the 

index of the superfield / is shifted by +1 (—1), i.e., fk = fk+i (fk = fk-i), an d it denotes the 
action of the direct (inverse) f-Toda transformation applied to the superfield /. The invariance 
of the N = 2 super-NLS hierarchy with respect to the transformation (|2.4j ) can be checked 
directly, and it means, in particular, that the right-hand sides of eqs. ( |2.1p are the solutions of 
the symmetry equationP] corresponding to the mapping (|2~~4 ) |l] . 



1 The sign ' means the derivative with respect to z. 

2 Let us recall that the symmetry equation for a given system can be obtained by differentiation of the system 
with respect to an arbitrary parameter. 
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3 The component form and boundary conditions 



The N = 2 superfield form (|2.4j) of the f-Toda mapping is not very suitable for actual calcula- 
tions and we rewrite it in terms of the superfield components defined as 

Vk = Dj k \ 1 ip k = Dj k \, u k = Df k \, ip k = Df k \, (3.1) 

where the index k 6 Z and | means the (8, 6) —>■ limit. In terms of these components, eq. 
( |2.4| ) becomes the infinite-dimensional chain of equations 

(\n(u k+1 v k )) ' = ip k+1 ip k+1 - ip k ip k} (3.2) 



-v k+1 ^ k+1 + v k ip k} (3.3) 
u k+i*P k +i ~ u k^ k i ( 3 - 4 ) 



Vk 

- (InUfc) " = u k+1 v k+1 - u k v k - ip k +i ' 4> k+l + 4> k ' ip k . (3.5) 

Each point of this chain contains four functions: two bosonic u k , v k and two fermionic tp k , ip k 
functions. In the bosonic limit (i.e., as ip k — > and ip k — > 0) it is equivalent to the usual Toda 
chain| ||. This is the reason why the chain ( |2.4|) is called the fermionic Toda or f-Toda chain 

E- 

For completeness, let us also present the component form of the global N = 2 supertrans- 
lations 

-/!_,- - 1 _- / 

v k -> v k - eip k - -eevk , ^ k -»• W k + + , 

1 1 
Uk-^Uk- eip k ' + 2 £eM fc '■> ^k + zu k - -ss4>k ' (3.6) 



under which eqs. (|3.2j )- (|3.5| ) remain invariant. Here, e and e are constant fermionic parameters. 

Equations (|3.2| )- (|3~5| ) also possess the inner discrete automorphism a m (m G Z) with the 
properties 

VmUkVm = V m~k, CTmVk^m = u m-k, CTm^k^m = ^m-ki ^m^k^m = 4>m-k, (3.7) 



which will be useful in what follows. Eqs. (|3.7|) can also be easily represented in terms of the 
superfields and superspace coordinates 

Vm/kCrn 1 — fm- k i a mfk a m = fm-k, cF m za^ = z, a m 6a^ = 6, o„fio~l = 9. (3.8) 



Let us discuss the boundary conditions which can supplement eqs. ( |3.2j )-(|3T 
We call the system ( p.2|) -( ^75|) with additional manifest N = 2 supersymmetric boundary 
conditions 

fo = & u = ip = 0, k > 0, (3.9) 

or 

J M = 0^v M = ^ M = Q i k<M, (3.10) 



3 In the bosonic limit, eq. (|3.2| ) admits the obvious solution u k +\ ~ — which, being substituted into eq. 
(^5|), produces the well-known representation for the usual Toda chain equations. 
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where M G Z + , the f-Toda chain interrupted from the left or from the right, respectively. The 
case where these conditions are satisfied simultaneously corresponds to the f-Toda chain with 
fixed ends. 

Application of the transformation <jm ( |3.7|) to the f-Toda chain interrupted from the left 
transforms it into the f-Toda chain interrupted from the right, and vice versa. Due to their 
relation, it is sufficient to consider only one of them; for concreteness, we discuss the f-Toda 
chain interrupted from the left in what follows. As concerns the f-Toda chain with fixed ends, 
it is invariant with respect to om-, or, in other words, it possesses the inner automorphism a M . 



4 The integrability of the f-Toda chain interrupted from 
the left 

In this section, we construct the recurrent procedure for solving of the f-Toda chain equations 
(|3.2| )- (|3.5| ) with the boundary conditions (|3.9[) . 



4.1 The conservation laws 

Our first goal is to resolve the f-Toda chain interrupted from the left. In other words, we assume 
that the boundary condition (3.9) is satisfied and Vq,i/j are arbitrary functions. The problem 



consists in expressing, using the f-Toda chain equations, all of the functions Uk, Vk, ?pk and ip k 
at k > 1 in terms of the functions Vo and ip , as well as a number of additional constants which 
can arise in integrating the equations. 

Let us start with a discussion of the conservation laws which are relevant to the problem 
under consideration. 

By a direct check, one can verify that the following set of conservation laws: 

Cfc-i = MfeWfc-i Wk-i + WkVk-i (4-1) 

Vk-l 

takes place at each value of k, where c& is an arbitrary constant. However, there is a simpler 



way to check this statement by rewriting (4.1) in the superfield form 



_ DDjfkh^Df^) _ DD(fkfk-iDfk) (A2) 
* 1 Dfk-i Df k 1 ' ) 



using eq. Q2.4|) and definitions (|3.1| ). The last equality of formula ( |4.2|) demonstrates that Ck is a 
constant (i.e., it is the integral of motion), because it is both the chiral and antichiral superfield, 
i.e. Dck = Dc k = 0. The action of the transformation a m ( |3.7| ) interchanges different integrals 
in accordance with the following law: 

O'mCfcO-" 1 = Cm-k-i. (4.3) 

Using the conservation law ( [4. 1[ ) , eq. ( |3.3j ) for the fermionic function ip k can identically be 
rewritten as _ 

(Vfc-l^fc) ' , ^fc-l^fc-lx I 1 / , (A .K 

( ) -«fc-i|Ufc-i - ( ) \Wh = -Ck-xWk-x- (4-4) 

Ufc-i ^fc-i 

From eq. ( f4.4|) , we observe that if we want to express ipk as a functional of the chain functions 
Ui, Vi, ipi and ^ defined at the previous (from the left) points of the chain (i.e., at i < k — 1), 
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it is necessary to solve a linear equation of the second order with the coefficient functions that 
also depend on the functions of the previous points of the chain. If we could find such an 
expression for ip k , all other independent functions ip k , u k and v k of the f-Toda chain could also 
be obtained in a recurrent way as functionals of the previous points of the chain. Indeed, using 
eqs. ( j3.4|) , (|3.5|) and conservation law (3~T), one can easily obtain the corresponding formulae 
for the functions i[) k , v k and u k . 

T I /^fc-l'w, T Ml , K-lV'k) ' T , -r 's 

Wk = {-Vk-l{ ) + Wfc-lV>fc_i)/(Cfc_i + 4>k-l - Wfc-i , 

Vk-l Vk-1 

v k = u fc _ 1 (-(ln + n fc _ii; fc _i + ipk^k ~ ^fc-iVfc-i)/ ( c fc-i + ^ kVk ~^ ^) fe _ 1 - ^ fe Vfc-i')> 

{^kVk-l) ' —r , —r /\ / , . 

M fc = (C fc _! + -0fc-l - Wfc-1 )M-1> ( 4 - 5 ) 

respectively. 

Thus, the crucial problem is the solution of eq. ( |4 . 4| ) for ^ by some recurrent procedure. 
This task is solved in the next subsection. 



4.2 The integrable factors 

The aim of this subsection is to construct the recurrent procedure for solving eq. ( |4.4| ) for an 
arbitrary chain point. 

Let us assume that the homogeneous part of the second order equation ( (4.4| ) possesses a 
bosonic integrable factor fi k . This means that after multiplication by ///-, the corresponding 
equation may be represented as a whole derivative. Following this line, let us rewrite eqs. ( |4.4| ) 
in the following equivalent form: 

/ Uk 2 (Vk-llpk \ I \ I r / ' \ / . r /^fc-lV'jfe-l \ I n / / lAR\ 

( ( ) ) -Uk-i{-( ) +A*fcWfc-i-( ) \}iPk = -Cfc-iWfc-i- (4.6) 

Equating the coefficient of the function ip k in the l.h.s. of eq. ( |4.6j ) to zero, we obtain the 
following equation: 

( ) = «fc-i - ( ) W (4.7) 

for the integrable factor fi k . At k — 1, it has an obvious solution /ii = 1 because, in this case, 
its r.h.s. becomes equal to zero due to the boundary conditions ( |3.9|) . 

If the integrable factors are restricted by eqs. (|4.7| ), eq. (|4.6Q for the function ip k becomes 



Vk-l Uk 

and one can easily integrate it, 



-Cfc_i/i fc Vfc-i (4-8) 



i>k = — — [-c fc _i / cte-^ / dzfi k ip k -i + Pk dz^—^- + a k ], (4.9) 
Vk-i J n k 2 J J fi k 2 

where a k and f3 k are fermionic constants of integration, and one can find the explicit expression 
for ip k if all other parts of formula (|4.9| ) are known. The appearance of these additional con- 
stants is a qualitative difference of the f-Toda chain in comparison with the integrable chains 
investigated before. 
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Concerning eqs. ( |4.7| ) for the factors [i kl they can be resolved in recurrent form with the 
following answer: 

/i fc+ i = w fe _i( ) ( ) ( ) , /ii = 1. (4.10) 

Indeed, this can be checked by direct substitution, taking into account the f-Toda chain equa- 
tions (|3.2| )-( pT5D . If the integrable factor ji k satisfies eq. ( |4.7|) , the factor Vk+i ( fll0|) satisfies the 
equation 

fVk + 1 \ I r /VV0fc\ / 1 (a i -i \ 

( ) =[«*-( — ) \vk+i- (4.ii) 

Vk V k 

This verification can be simplified if one uses the following convenient identity: 

(— (u fc _i( ) ) ) =[ix fc -( ) j u fc _i( ) -c fc _i( — ) , (4.12) 

which can easily be checked. 

Let us stress that the integrable factor Vk+i ( 4.10| ) is the functional of the chain functions 
defined at the previous (from the left) points of the chain, thus, the same important property 
is also satisfied for the function ip k fl4.9|) . 

To consider the case of the f-Toda chain interrupted from the right with the boundary 
conditions ( |3.10| ), one can simply apply the transformation <jm ( |3.7|) to formulae (|Q|)-( fLT2"| ). 
Without additional comments, let us present expressions for ip k and their integrable factors ~p k 

i>k = -^-[-Cfc / dz^^ / dzJL k i) k+l + p k dz^P^- + a k ], (4.13) 
Uk+i J Vk J J Vk 

— I Vk \ I , Vk (4>k + l, , ( Uk+X'4 ! k\ I — A I A -.AS 

V k -i = u k+1 { ) + ( ) (— ) , /i M _! = l, (4.14) 

U k +l C k U k+ \ fi k 

where a k and j3 k are fermionic constants, which will be useful in what follows. 
Thus, the following proposition summarizes this section. 

The solution of the f-Toda chain interrupted from the left is given by recurrent relations 
(0, gJfy and (O. 



5 The integrability of the f-Toda chain with fixed ends 

In this section, we construct the recurrent procedure for solving the f-Toda chain with the 
boundary conditions ( |3.9|) , (|3.10|) . 

5.1 New form of the f-Toda chain 

Let us introduce the new basis q k , r k , £ k and £ fc in the space of the f-Toda chain functions u k , 
v k} ipk and ip k , defined by the following invertible transformation: 

k k 

Qk = n( M ^-i)' r o = ^o, r k = v k If ( u sV s -i), 

s=l s=l 

k k 
?0=^0> Ck = TpkTKVaVs-l), 6 = t/jk/Y[(u 8 V 8 -i), (5.1) 

8=1 8=1 
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Uk 



Tk-l 



v = r , v k 



r k -j T T £fc 

— , = £0, ^k = —■ 

qk qk 



i>k = Cfc9fc> 



(5.2) 



where k = 1,2, ...,M. In this case, the boundary conditions ( ^.9| ) and ( |3.1C| ) are transformed 
into the following conditions: 

go = £o = 0, (5.3) 



M 



0. 



(5.4) 



respectively. In the new basis, the f-Toda chain equations (|3.2j) - (|3Tg) have the following form: 



On 7,)' Os,- (lng M ) ' = 0, 



(5.5) 



\Jm~\Km ) = rM-KM-i, -( ) = — , 

r r 



r A + O is 



(5.6) 
(5.7) 



r .i-i 



(5.8) 



-(lnr )"- ri * ,T n^w^+i 



r 



(5.9) 



Hereafter, the index j lies in the following range 1 < j '• < M — 1. 



The substitution of the transformations (|5.1|) , (|5.2|) into (|3.7|) at m = M gives the inner 
automorphism cr M , 



qu 
qu-k 



qu 



^M-k-l 



OhlKkO M 



M-k 

qu 



^M^M = QM^M-k, (5-10) 



foreqs. (PHPD- 

It is interesting to note that the functions qk, Cat an d £ are completely decoupled from eqs. 
(|5.7|) , (|5.8|) and (|5.9|) , which form a closed set of equations for the functions r , tj ^ and 



Moreover, taking into account the corollary from eq. (|5.5| ) that g^/ is an arbitrary constant, 
one can conclude that the transformations ( |5.10| ) are also closed for them. However, this is not 
the case for the N = 2 supersymmetry transformations (|3.6|) , which, for eqs. (|5.5|) - (|5.9|) , have 
the following infinitesimal form: 



Sqk 



c - i j_ Vkdqk Ik^qk Xi 

+ > 5 ik = £r k + — — , 



r k -i 



qk 



r k -i 



£k5qk 
qk 



(5.11) 



where eqs. (|5.5| ) have been used. 

We call eqs. (|5.7|)- fl5T§D the restricted f-Toda chain. If some of their explicit solutions are 
known, one can also easily obtain the solutions of eqs. Q5.5|) -( [5l)D for qj, £ M and £ , 



fat 



£, + £ i+1 £/), gi = cft/(l 



(r 6) V 



r 



to + teo'), (5-12) 
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£m = - / dz—^— [ dzr M -i^M-i + Pm [ dz— h a M , 

J r M -x J J Tm-x 

£ = - J dzr J dz-^- + Po J dzr + a , (5.13) 

where aM,«o an d Pm, Pq are the fermionic constants of integration and we have used the 
conservation laws (|4.1| ) expressed in terms of the new functions qk, r^, and £ fc ( |5.1|) , ( |5.2|) . 
Let us introduce the new bosonic function Xj 

r = exp(x ), Tj = exp(xj), exp(— x_i) = exp(x M ) = 0, (5-14) 

which is usually used in the case of the bosonic Toda chain and is more suitable for the system 
under consideration. In this case, the restricted f-Toda chain equations ( |5.7|) - (|5T9|) and their 
inner automorphism <7m become 

" + Xj-i ' Cj ' - exp(xj - xj-!) = (5.15) 



f j " ~ x i ' €j ' ~ exp(xj - x^) = (5.16) 
x " + exp(a;i - x ) - £i ' f x = 0, 

Xj " + exp(x j+ i - Xj) - exp(x i - Xj- X ) - £j +1 ' - ^ ' = 0, (5.17) 

(T mXq& M = —XM-Xi OM x j&M = ~ x M-j-l, <J M^j cr AI = ^M-ji n Ws /' A/' = ^M-j, (5.18) 

respectively. In what follows, we concentrate on the analysis of these equations and prove their 
integrability. 

5.2 The general assertion 

In this subsection, we briefly discuss the general problem of integrating a system of ordinary 
differential equations containing both unknown fermionic and bosonic functions. 

Althought this is not crucial for the general discussion, for defmiteness, we assume that the 
system consists of the second order differential equations for M bosonic and 2(M — 1) fermionic 
independent functions, as takes place for eqs. ( |5.15| )-( [5.17| ). In this case, its general exact 



solution must include 2M bosonic and 4(M — 1) fermionic independent constants, of course, if 
the system is integrable. 

Let us analyze the qualitative structure of this exact solution in more detail. 

Evidently, the exact solution is a polynomial with respect to fermionic constants for any 
function involved in the system. It is well known that the bosonic (fermionic) functions are 
even (odd) polynomials and any fermionic constant enters a monomial representing the product 
of fermionic components only once[[ 

It is instructive to treat this exact solution as the result of some iteration procedure or as 
a result of calculations in the framework of some perturbation theory, where the role of per- 
turbation parameters is played by the fermionic constants, and different orders of perturbation 

4 Let us recall that each polynomial is the sum of monomials composed of the products of different fermionic 
components and, by definition, its degree is the number of components of its maximum monomial. 
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calculations are in one-to-one correspondence with different degrees of the monomials, i.e., with 
the different numbers of fermionic components composing them. 

Now we will analyze the qualitative structure of different orders of perturbation calculations 
in the framework of such a perturbation theory. 

The zero-order approximation to the exact solution corresponds to the general solution of 
a pure bosonic system, which can be derived by taking the bosonic limit of the initial system. 
All of the bosonic integration constants parametrising the exact solution of the initial system 
must arise in this order of perturbative theory. 

The first order approximation corresponds to the general solution of a linear fermionic 
system, which can be derived by substituting the zero-order approximation for bosonic functions 
into the initial system linearized with respect to fermionic functions. Similarly to the case of 
zero-order approximation, all fermionic integration constants parametrising the exact solution 
of the initial system must arise in this order of perturbative theory. 

The second-order correction corresponds to the solution of a linear bosonic system with the 
following structure: its homogeneous part is the symmetry equation 2 of a zero-order bosonic 
system and its inhomogeneous part is what remains after the substitution of the zero and first 
orders of perturbative calculations into the initial system and rejection of the monomials of all 
degrees but monomials of the second degree. 

The homogeneous part of a fermionic system, corresponding to the third order-correction 
to the exact solution, coincides with the fermionic system of the first order- approximation and 
its inhomogeneous part is given by perturbative decomposition of the initial system using the 
functions of the first and second orders of the perturbative calculations, and so on. 

To close this general discussion, let us stress that the homogeneous part of a system corre- 
sponding to any even (odd) order- correct ion to the exact solution coincides with the system of 
the second (first) order- approximation and neither new bosonic nor fermionic independent pa- 
rameters appear in any order starting with the second order of the perturbation calculations^]. 
We would also like to emphasize that by construction, such a perturbation theory is convergent 
and it gives the exact solution: the perturbation series is interrupted due to the fermionic na- 
ture of the perturbative parameters because a limited number of monomials can be composed 
using a finite set of fermionic constants (i.e., the solutions of the first-order approximation). 

Now let us formulate the assertion with respect to the initial system: 

if the equations of the zero and first orders of the perturbation theory are exactly integrable, 
the initial system is also integrable, at least in the quadratures. 
The proof is given in a few words. 

If the equations of the zero (first) order of the perturbative theory are exactly integrable, 
the corresponding symmetry equation is also exactly integrable. Indeed, one can derive its 2M 
(4(M — 1)) linear-independent solutions by taking the derivatives of the general solution of the 
zero (first) order equations with respect to its 2M (4(M — 1)) bosonic (fermionic) indepen- 
dent constants. As follows from the previous discussion, the problem of resolving equations 
corresponding to any other order of the perturbative calculations reduces to the problem of 
resolving the inhomogeneous symmetry equation, but the last problem is an exactly solvable 
one. Indeed, if the general solution of some linear homogeneous system is known, like in the 
case under consideration, then, by applying the well-known method of varying the arbitrary 
constants, one can algorithmically construct the solution of the corresponding inhomogeneous 

5 Clearly, having arisen, they can always be eliminated by redefining the parameters of the zero and first 
orders of the perturbation calculations. 
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equation, at least, in the quadratures. Thus, all the orders of the above-discussed perturbation 
calculations can be resolved in explicit form. 



5.3 Solution of the restricted f-Toda chain 



The purpose of this subsection is to apply the above-developed regular algorithm of integrating 
a system of ordinary differential equations containing both unknown fermionic and bosonic 
functions to the restricted f-Toda chain equations ( |5.15|) -( |5.17|) . 

Following the line of section 5.2, let us represent the functions xq, Xj, £j and £ • as a pertur- 
bation series, 



2(M-1) 
1=0 



2(M-1) 

E 

1=0 



£3 



2(Af-l) 

i=i 



2(M-1) 



(5.19) 



where the perturbation parameters coincide with the 4(M — 1) fermionic constants of the 
first-order approximation to the general solution. Here, the functions x^ and x^ (£ J - 2 ' +1 ' ) 

and ) are the 2/ (2/ + 1) order corrections to the zero (first) order approximation. The 

fermionic character of the decomposition parameters guarantees that this series is interrupted, 
starting with the 4(M — 1) + 1 order. As result, such a perturbation theory is convergent and 
gives an exact result. 

Substituting decompositions ( |5.19| ) into eqs. ( |5.15| )-( |5.17| ), extracting terms of the same 
order, and equating their algebraic sum to zero independently for different orders, one can 
obtain the complete set of perturbative equations corresponding to the presented perturbation 
theory. 

Here, we demonstrate that the conditions of the general assertion in the previous subsection 
is satisfied for the restricted f-Toda chain, i.e., it is exactly integrable. 

The zero order bosonic system coincides with the usual one- dimensional Toda chain || 



(°) " _L amtaM J°) 



0. 



X 



+ exp(x [ f_l_ 1 — Xj 1 ') — exp(x^ 0) — a^i) = 0, 



,(o) 



(5.20) 



which is exactly integrable. Its general solution is well known and can be represented in the 
following form ||: 



,(o) 



( (0)\ 

exp(xo ) 



Al 



Oj exp (bix), exp (xj 



(0), 



' Deb* 



(5.21) 



1=1 ^^Vj 

where a>i and bi are arbitrary constants, and Detj is the jth principal minor of the matrix 



( r o 



(0) 



r 



Jo) / 



(0) 




r, 



(0) // 




,(0) // (0) ,// 







r 



(0) // „(0) /// 













,(0) 



V 



...\ 



/ 



(5.22) 



The first order fermionic systems have the following form: 



+ x 



(0) / / 
3-1 



^'-^exp^-o^O + e 



(o) 



.(0) 



r(i) 



(5.23) 
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(!) // _(0) / 7(1) / 7(1) _/,JO) _(0) x , -(1) 



-X) 



^r , -6'«p(4 ; - a? i-i)+e+i=° 



(5.24) 



These equations are also exactly integrable and using the results of section 4.2, one can easily 
generate their general solutions. Thus, substituting transformations ( |5.2|) , (|5.14|) , as well as 
solutions ( |5.12[ ), into relations (|4.9|), ( |4.1U| ), ( |4.13| ), and ( (4.14[ ), and rejecting terms nonlinear 
with respect to the fermionic fields, we get the following recurrent formulae^: 



exp(xfl 1 



(0)' 



+ a j ], 



// 0) -/y (0) 



T (0) / ..(0) 
•^j— 1 o ' 



/y (0) 



(5.25) 



^=exp(4^f[- 



dz 



exp(— x 



(Oh 



dzfif^ + /3j / dz 



exp(— Xj 



(0)- 



+ a 



(5.26) 



for the exact solutions of eqs. (|5.23| )- (|5.24j) . Thus, following the general assertion of section 
5.2, we conclude that the restricted f-Toda chain is also integrable. 

Now we show how to construct the solutions of the perturbation equations for all other 
orders of the perturbation theory with respect to the 4(M — 1) fermionic constants atj, j3j, aj, 
and of the first order solutions ( |5.25| ), ( |5.26| ). 

For the perturbation equations of the 21 (I > 1) and 21 — 1 (/ > 2) orders, we have 



x 



(21) „ 


(20 „ 



x 



(21) 
1 

.(20 



— x 21 ^) exp(xi 



(0) 



(0)> 



X, 



+ ( x j+i ~ x f lS> ) ex ${ x Y+i — x j' > ) ~ i x Y l> — exp(x] u '' — xfli) = X 
J2i-i) „ (o) / .(2i-i) , A2\-\) ( (o) (o) v .(a- 1 ) _ 



,(0) 



,(0) 



(20 

5 

(20 



.(20 ■ 



„(°) 



.(0) 



'(20 



(5.27) 
(5.28) 



-(2i-i) ,/ 
4, - x • 



^'ef^'-ef^exp^f 



( 2 „„J v(») fwW- 1 ) 



X 



(0) • 

i-i. 



(21-1) — (21-1) 



(5.29) 



and H - 2Z ^) are defined by the following 



respectively. Here, the functions Xq and Xj (S^ 

' — - exp(xi - x ) + £i ' £i + (1 + x 1 — Xq — Xi 0) + Xq 0) ) exp(xi 0) - Xq) 



relations: 

2(M-1) 



x = 



X, 



E ^ 2 

2=1 
2(M-1) 

E x f l) 
i=i 



exp(x j+ i - xj) + exp(xj - x^i) + ' £ j+x + ^ ^ ' 



+ (1 + x j+i x j x j+i x j ) e -^p( x j+i • i 
(1 ~\~ Xj Xj—\ Xj ~\~ Xj _. j ) exp(x^ ^) 

2(M-1) 

_ j2 Ef l ~ 1] = (xf\ - Xj_i) ' + £j(exp(x 

1=2 
2(M-1) 



„(0h 



„(0) 



(o)^ 



— exp(x!j 0) — XjHi)) 



(o) 



• (a;.- 



' ' + ^(exp(x j - - exp(x}° } - x^)) 



1=2 



3 Here, we have rescaled the fermionic constants of expressions (4.9) and ( 4.13 ). 



(5.30) 
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respectively, and they are the sum of the monomials of 21 (21 — 1) degreeQ Let us stress 
that the inhomogeneous parts of perturbative equations ( |5.27|) - (|5.29|) , corresponding to a given 
order n, depend only on the functions XQ m \ xf m \ ^ 2l+1 ^ and of the previous orders of 

perturbative calculations (i.e., at m < (n — l)/2 and / < (n — 2)/2). Thus, we have a consistent 
perturbative theory. 

If we set the inhomogeneous parts of perturbative equations ( |5.27[ ) equal to zero for any 
even order of the perturbative theory, they coincide with the symmetry equations 2 for the 
Toda chain ( 5.20|) . Thus, in this case, their 2M linear-independent solutions yfc at each point 
k (k — 0, 1, M — 1) of the chain are given by 



dx 



(0) 



dAf 



(5.31) 



where the capital Greek letter-indices run over the range A, $ = 0, 1, ...,2M — 1, and A\ = 
{a\, clm-i bi, bu } (see eq. (|5.21| )). To use the formalism of varied constants for the solution 
of inhomogeneous equations (|5.27 ), it is necessary to solve the linear system of the first order 
differential equations 



2M-1 2M-1 

E vtcf\z)' = ^ E yt'cf\z)' = x^ 

A=0 A=0 



(5.32) 



with respect to (z), where ^'(z) are 2M bosonic parameter- functions. Then, the solution 
of eqs. ( 5.27 ) for the function has the following form: 



,(20, 



r (20 



2M-1 

E vtcf\z). 

A=0 



(5.33) 



Let us introduce the 2M x 2M matrix V§ defined as 



v 



yl 



(5.34) 



Then the solution of the system ( |5.32| ) can be represented in the form 



,(20/ 



M-1 



fdz E (P- 1 )^, 

J m=0 



(5.35) 



where V 1 is the inverse matrix for the matrix V (V 1 V = VV 1 = I) and expression (|5.33| ) 



for x^ becomes 



x 



(2/) 



2M-1 



M-1 



A=0 m=0 



(5.36) 



Taking the sum over all orders of the perturbative theory and using relations (|5.19|) and ( |5.30|) , 
we have the following exact expression for Xk'- 



(o) 

•"k ^k 



2M-1 



M-1 



+ E Pu +k jdzY, {P~ l )tx m . 

A=0 m=0 



(5.37) 



7 It is implied that the perturbative decompositions (5.19) of the functions xq, xj, £j and ^ - must be substi- 
tuted into the r.h.s. of relations ( 5.30| ). 
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One can invert the above-developed perturbation scheme and consider the relation ( |5.37 ) as 
the equation for Xk, and take it as a starting point. Then its iteration (|5.19| ), together with 
the iteration of the corresponding equations for the fermionic functions £j and £ • (see below), 
is interrupted starting with the 4(M — 1) + 1 order and gives the exact solution for X}.. 

The odd-order perturbative equations ( |5.28|) , ( |5.29| ) can be integrated in the same way as 
described for the case of the even-order equations ( |5.27|) . Without going into detail, let us 

present only the expressions for the functions ^ 2 ' -1 ^ and ^ 1 , 

2(M-1) rj e (l) 2(M-1) nT^ 1 ) 

«f" = £ S-cT"^), ??'-"= £ S-JS-"^), (5-38) 
n=i 0,1 n n=i 0,1 « 

where T n = (a x , a A /-i, A, /?a/-i), F n = (a%, ...,a M -i, Pi, P M -i) ( fi = 1,2, 2(M - 
1)); Cq' ^(z) and Cq 1 (z) are 4(M — 1) fermionic parameter-functions that are solutions of 
the following linear system: 



2(M-1) Qt(l) 2(M-1) o^(l) 

0-! ^ W " ' ^ 



V ^-^'^(z) '-n V 



2(M-1) otC 1 ) 2(M-1) ^(1) 

£ §-<r ' = 0, E (§") ' ' = 3f -". (5.39) 

o=i 0,1 o n=i 0,1 o 

To close this subsection, we would like briefly discuss a slightly modified version of the 
perturbative theory for the odd-order perturbative equations ( |5.28|) and (|5.29| ). 

To do this, let us represent the differential equations ( |5.28| ) and ( p.29| ) in integral form. Thus, 
multiplying them by the integrable factors (if 1 and Ji f , respectively, using the equations 

(eM-xf\) Mf ' ) ' = eM-xf-*) («p(*f ) Jlf ')' = exp(zf x ) /if (5.40) 

for the factors, taking the sum of all perturbative orders, and using relations ( |5.30|) , eqs. (|5.28|) - 
(|5.29 ) can be identically rewritten in the form 

( (°) A 

/ / (o) \ (o) 2 ,exp\Xj-i) t ww (o)/ t , ^ n 
{exp{-x)Y) fi) ( ^ Q ) = fi) + ^), 

w 

w 

which can be easily integrated. As a result, we have the desirable integral form of eqs. 
(|5.29 ) given by 

c ( (°) ^ (°)r I a ex P(4-i) [j (o)/ t , ~ n , fl /" . expjxfp 
& = exp(-x}i 1 )/i} ; [ / dz -4— / dzfi) + a,) + Pj / ^ ^p— + a,], 

ny Yj 



-P ( (o)s_(o) r / , exp(-xj 0) ) r _ (0) r exp(-xj 0) ) _ 

f j = exp(a;} >} [_/ tfo — _ (0)2 _/ ^ ("fi+i + -i) + 0j / dz — _ (0) 2 + a iJ- ( 5 - 42 ) 

w w' 



Simple inspection of eqs. Q5.42 ) shows that their solutions can be consistently obtained by iter- 
ating with respect to the fermionic constants of integration a,, Pj, oij and Pj in the framework 
of the above-discussed perturbative scheme. 

Thus, all the orders of the perturbation calculations can be resolved in explicit form. 
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5.4 Example: the M = 2 case 

To illustrate the general formulae of the previous subsection, here we consider the simplest 
example of the restricted f-Toda chain with fixed ends at M — 2. 

In this case, eqs. ( |5.15| )-( |5.17| ), boundary conditions (|5.3|), ( |5.4| ), (|5.14j) , as well as pertur- 
bative decompositions (|5.19|) , have the following form: 

£i " + x ' £i ' - & exp(xi - x ) = 0, ?]. " - a* ' £ t ' - Ci exp(x x - x ) = 0, 

x " + exp(x 1 - x ) - £i ' £i = 0, xi " - exp(a;i - x ) - £i £i ' = 0; (5.43) 

£o = £ 2 = exp(-s_i) = exp(x 2 ) = 0; (5.44) 

* = 4 0) +4 2) +4 4) , x 1= x^ + 4 2) + xS 4 \ ^=d 1} +d 3) , f 1 =?f ) +ff ) , (5-45) 
respectively. 

According to formulae ( |5.21|) , the zero-order functions Xq and xf^ are given by 
exp(xo°' ) ) = aiexp(biz) + a 2 exp(o2<2), exp(x^) = — 0102(61 — O2) 2 exp(z(6i + 62) — x< o^)- (5.46) 

Substituting them into eqs. (|5.25|) - ([5.26|) , taking into account eqs. (|5.44|) and integrating the 
obtained expressions, we obtain the following results: 

= foxf ' + a x exp(-xi 0) ), = P.x^ 1 + a x exp(4 0) ) (5.47) 

for the first-order functions and ^ . 
According to formulae ( |5.30|) , we have 

v (2) _ A\) / 7(1) y (2) _ Ai) / 
A o — Si Si ) A i — Si Si ) 
=(3) = _ d D , 42) , + g )[x?) _ x ( 2)) exp(x ,„) _ ^ 

sf>=i ( ,"'4 2 "+?i I, (4 2) -4 2 »)exp(xr-4° ) ), 

A-r = -i(4 2> - 4 2, ) 2 exp(4° , - 4°') + d 1 ' ' If + d 3) ' li 1 ', 

Af> = i(4 2 » - 4 2 ') 2 exp(4 0) - 4°>) + «P f f ' ' + «?> fS" ' (5.48) 

for the inhomogeneous parts of perturbative equations (|5.27|) -( |5.29|) , corresponding to the 
higher order corrections. Iterating eqs. (|5.42p and (|5.37| ) in consecutive order, step by step, we 
obtain the following expressions: 

4 2) = - n^i^i ~~ "iai)4 0) ' - ai/3iexp(-xj 0) ), 

X W = -{p{j} x - aiai)x[, 0) ' + Aai exp(x^) (5.49) 
for the second-order functions, 

d 3) = ~ai/3i^(exp(-4 0) )) ' + «i/?i«i(^ exp(xS 0) - xf) - bfo), 

ti ] = -«i«i /3 1 (-exp(xS° ) - x[, 0) ) - 6i6 2 ) - /?i(exp(xS 0) )) ' (5.50) 

for the third order functions, and, at last, 

4 4) = -xf ] = -ocxPioix exp(xS 0) - xi 0) ) - &i& 2 ) (5.51) 
for the fourth order functions. 
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6 Conclusion 



In the present Letter, we proved the integrability of the f-Toda chain with fixed ends or inter- 
rupted from the left (right) and proposed an algorithmic method for constructing its explicit 
solutions. Many interesting questions arise in this connection. What is the group-theoretical 
foundation of this result and its connection with the representation theory of supergroups (su- 
peralgebras)? Is it possible to construct a superintegrable two-dimensional generalization of 
the f-Toda chain? Is it possible to represent the f-Toda chain with fixed ends in the Lax-pair or 
Hamiltonian forms? This last question is a very important in connection with the problem of 
its quantization. In the case of the usual Toda chain, all these questions have the answers and 
the authors hope to find their solutions for the case of the f-Toda chain in future publications. 
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